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ABSTRACT 

We present a new analysis of the Aquarius simulations done in combination with a 
semi-analytic galaxy formation model. Our goal is to establish whether the subha- 
los present in ACDM simulations of Milky Way-like systems could host the dwarf 
spheroidal (dSph) satellites of our Galaxy. Our analysis shows that, contrary to what 
has been assumed in most previous work, the mass profiles of subhalos are generally 
not well fit by NFW models but that Einasto profiles are preferred. We find that for 
shape parameters a = 0.2 — 0.5 and v max = 10 — 30 km/s there is very good corre- 
spondance with the observational constraints obtained for the nine brightest dSph of 
the Milky Way. Furthermore, the internal dynamics of these systems, as well as the 
number of objects of a given circular velocity are also matched if the total mass of the 
Milky Way is ~ 8 x 10 11 M , a value that is in agreement with many recent deter- 
minations. Our simulations show important scatter in the number of bright satellites, 
even when the Aquarius Milky Way-like hosts are scaled to a common mass, and we 
find no evidence for a missing population of massive subhalos in the Galaxy. This 
conclusion is also supported when we examine the dynamics of the satellites of M31. 
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1 INTRODUCTION 

Despite the great success of the ACDM concordance cosmo- 
logical model on large scales, on the scales of galaxies and 
below the theory is often defied. Some of the issues on small 
scales have been consistently explained within the theory 
itself with the inclusion of physical processes that mostly 
affect baryons. This is the case for the "missing satellite 
problem" (Kl ypin et al.|[T999[ ), namely the overabundance 
of satellites in dark matter only simulations compared to 
the observed number of luminous objects around the Milky 
Way and other nearby galaxies. It is now widely accepted 
that the shallow potential wells of small dark matter ha- 
los must be strongly affected for example, by reionization, 
feedback and that star formation is highly inefficient in such 
systems ( Somerville|2002 ' I. 

Recently, |Boylan-Kolchin, Bullock fc ICa plinghat 
( 2011a|b I have argued that the dark matter satellites (sub- 
halos herafter) predicted by ACDM are persistently too 
dense to host the observed population of dwarf spheroidal 



E-mail: cavera@astro.rug.nl 



galaxies (dSph) in the Milky Way if these are embedded in 
halos following |Navarro, Frenk fc White| ( |1996| |1997[ here- 
after NFW) profiles. A few studies have been presented in 
the literature that address this conundrum. ILovell et al.l 
( |2011[ ) showed that in warm dark matter cosmological simu- 
lations of Milky Way-like halos, the circular velocity curves 
of subhalos are consistent with those derived by | Wolf et al.| 
(|2010[> for the Milky Way sate llites. Following a similar line, 
|Vogelsberger, Zavala fc Loeb (20121 carried out simulations 
of self-interacting dark matter and showed that the most 
massive subhalos develop cores, what could partially solve 
the problem. On the other hand, di Cintio et al. (20111 



pointed out that by including baryons in the cold dark- 
matter context the problem becomes more severe probably 
due to the additional adiabatic contraction experienced by 
the dark matter subhalos hosting gas. 

There are two assumptions implicit in the models which 
may lead to biased answers if overlooked. These concern (i) 
the actual mass of the Milky Way and, (ii) the density pro- 
files followed by dark matter satellites assembled in ACDM. 
The first issue has been addressed with a plethora of meth- 
ods leading to measures that, usually, are consistent with 
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a total mass of 0.7 - 2.0 x 10 12 M Q (|Wilkinson fc Evans 



1999} ISakamoto, Chiba fc Beers|2003| |Battaglia et al.12005 
20061 ISmith et al||2007l|Li fc Wmte||2008| |Xue et al.||200S 
Kallivayalil et al.||2009||Guo et al.||2010| |Watkins, Evans 



An 20101. The measurements suffer from uncertainties in 



the modelling as well as limitations in the kinematics of the 
tracers used. Therefore comparisons to simulations of Milky 
Way dark matter halos should take into account this uncer- 
tainty. 

On the second issue, namely the density profile of sub- 
halos, significant progress has been made, especially in re- 



cent years. Already Stoehr et al. ( 2002 1 found that the circu- 



lar velocity curves of subhalos in cosmological N-body simu- 
lations are more narrowly peaked (in a log-log plot) than the 
widely used NFW models. The outstanding numerical res- 
olution achieved in the latest of such cosmological iV-body 
simulations has enabled a closer examination of the shape 
of the density profile down to the innermost few parsecs of 
dark matter halos ( Springel et al.|2 008 Madau, Diemand & 



Kuhlcn 2008 1. Such studies have shown that Einasto models 
provide better matches to the density profiles found in the 



simulations than the NFW form( |Navarro et al.||2010 Reed, 
Koushiap pas fc Gao||2011 1 



In this paper, we will reconsider both these issues and 
establish how much they affect the conclusions drawn by 



Boylan-Kolchin, Bullock & Kaplinghat ( 2011b I. Like these 



authors we use the Aquarius halos, but we supplement the 
dynamical information provided by the simulations with a 



semi-analytic galaxy formation model (Starkenburg et al 



in prep I. One of the advantages of this approach is that it 



enables us to directly compare objects in the simulation with 
those observed. We introduce some relevant features of the 
simulations in Section [2] while in Section [3] we present in 
detail the results of our analysis. We draw our conclusions 
in Section |4] 



2 NUMERICAL PRELIMINARIES 

We use the simulations of the Aquarius project, six Milky 
Way-sized dark matter halos assembled in a background cos- 
mology consistent with the constraints yielded by WMAP-1. 
Each halo (labeled from A to F) was simulated at different 
resolutions, starting from a particle mass m p — 3.143 x 10 6 
Mq for the lowest and m p = 1.712 x 10 3 M Q for the highest 
resolution. In what follows we focus on the level 2 which is 
the highest level at which all the Aquarius halos were simu- 
lated (For more details see Springel et al.(2~008 l. 

Subhalos in these simulations are identified as bound 
overdensities with SUBFIND (Springel et al.||2001|). For each 



subhalo we compute the circular velocity profile as t> 2 (r) = 
Gm(r)/r, where m(r) is the mass enclosed within the spher- 
ical radius r. The maximum circular velocity « max is defined 
as the peak of the circular velocity curve, and is reached at 
position r max , i.e. w c (r max ) = u max . 

Numerical convergence is established by looking at 



the convergence radius as defined by Power et al. (20031 



Navarro et al. (2010) showed that the roots of the equation 



/200 n(r) 
8 In n(r) 



p(r) 



-1/2 



(1) 



correspond to different degrees of convergence depending 



on the value of the parameter k. In this equation n(r) is 
the number of particles enclosed within the radius r and 
p(r)/p c is the spherical density at this position in units 
of the critical value. Comparing the various resolutions of 



the Aquarius simulations, Navarro et al. (2010) showed that 



k = (7.0, 1.0, 0.4) correspond to deviations in the circular ve- 
locity profile of about (2.5%, 10%, 15%) respectively. Here, 
we use r conv («; = 0.4) for each subhalo. We also define the 
tidal radius rydal °f a subhalo as that which encompasses 
95% of the bound particles. The results from our definition 
show a good agreement with the output from more sophis- 



ticated expressions for the tidal radius (Tormcn, Diaferio fc 



Syer||1998 1 



To make more direct comparisons to the satellite pop- 
ulation of the Milky Way, we have also run a semi-analytic 
model of galaxy formation for all the Aquarius halos. This 



model is based on that originally de veloped by [K auffrnann 
et al.l (119991); ISpringel et all rt2001h; IDe Lucia, Kauffmann 



fc White| ( |2004[ ); |Croton et al| ( |2006| ); |De Lucia fc Blaizot 
(20071 and later modified to describe more accurately pro- 



cesses on the scales of dwarf galaxies ( Li, De Lucia & Helmi 



20101. The implementation used here also includes recipes 



for stellar stripping and tidal disruption. The resulting satel- 
lite luminosity function agrees well with that of the Milky 



Way as reported by Koposov et al. (20081 (see Section 3.31. 
Also the internal properties of the satellites, such as scal- 
ing relations, metallicities and star formation histories are 
in good agreement with those observed (for more details see 
|Starkenburg et al.|in prep] ). 



3 RESULTS 

3.1 About the density profiles 

It was recently reported that mass profiles of ACDM halos 



deviate from the NFW functional form. Navarro et al. ( 2010 1 



showed that a parametric model with a density profile with 
logarithmic slope described by a power-law (Einasto profile) 
provides better fits for objects of virial mass ~ 10 12 Mq. 
The power index a adds another free parameter, therefore 
the fits are expected to improve. Nevertheless, it was shown 
by |Springel et al.| ( |2008"| that even after fixing a = 0.16 the 



Einasto profile still yields much better results. The nature 
of the shape parameter a has been recently investigated for 
isolated objects with masses in excess of 5 x 10 12 M Q and 
the results suggest a deep connection with the pseudo-phase- 
space density distribut ion (jLudlow et al.|2011[ ). For smaller 
objects Springel et al. ( 2008 1 also showed that the mass pro- 



files follow much closer the Einasto than the NFW model. 
For both models, the mass enclosed within the spherical ra- 
dius r can be written as 



m(r) = 47rr_ 2 P-2(?(r/r_2), 



(2a) 



where r_2 is the radius at which the logarithmic slope of 
the density profile reaches the isothermal value and p—i is 
the density at that position. The details of each model are 
inherited by the function g, which takes the form 



Snfw(ie) = 41n(l + x) 



Ax 

T+x' 



(2b) 
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Figure 1. Spherically averaged circular velocity profiles v^(r) 

,,2 



semi-analytic model. Velocities have been scaled to v^ 2 = 4nvGp— 2r^_2- As already reported in Stoehr et al.|||2002k the velocity profil 



Gm(r)/r for the subhalos that are predicted to host stars by our 
2 



of subhalos tend to be more narrowly peaked than in the NFW form. The sample of subhalos has been grouped according to the best 
fit value of a, and plotted with different colors. An Einasto profile with the average value of a for each bin is overplotted. The residuals 
from the best-Einasto (NFW) fits are shown in the middle (bottom) panel, and in general are consistent with zero for the Einasto profile 
and exhibit systematic deviations from zero for the NFW case. In the column (a) = 0.24 both models yield similar results, which is 
naturally expected since the NFW equivalent is reached with a = 0.22. The systematic change of r conv /r_2 with a is a consequence of 
setting the convergence parameter k to a fixed value. 



for the NFW profile, and 

, , 1 /31na + 2-ln8\ f 3 2x\ , n , 

ffEinasto(a;) = - exp 7 — , (2c) 

a \ a J \a a J 

for the Einasto model. Here 7(0, x) is the lower incomplete 
gamma function. Although intrinsically different, these pro- 
files resemble each other for a « 0.22 in 0.01 ^ r/r-2 ^ 100. 
That means that objects that have a shape parameter close 
to this value are well fitted by either model. Fig.[T]shows the 
spherically averaged circular velocity profiles n 2 = Gm(r)/r 
for all the subhalos that are predicted to host stars accord- 
ing to our semi-analytical model. In total we calculate 20 
bins in the region r conv sC r ^ 0.9r t idai, we use this upper 
cutoff to ensure that our fits are not driven by tidal effects. 
All objects have at least 200 particles, but generally signifi- 
cantly more than 1000. For each of the plotted subhalos we 
calculate the merit function 

iv blna 

E = V (lnv 2 c (n) - ln^) 2 , (3) 

t — 1 

and minimize it against the free parameters of each model. 
We have deliberately chosen to use the cumulative mass 
instead of the differential profile since it is less sensitive 
to the shot-noise of each bin, as a consequence we can 
go to low number of particles, whenever the restriction 
n(0.9r t idai) — n(r colw ) > 200 is met. 

Fig. [l] shows the results of our fitting procedure. Here 
the subhalos have been binned according to their best fit a 
value. The three different columns show the subhalos that 



fall into each cv-bin, the average a within each bin is quoted 
in the top-left panel. Each curve in the top row has been 
conveniently normalized to the characteristic velocity v 2 2 = 
4nGp-2r_2 an d the characteristic radius r_2. We have also 
overplotted the predicted Einasto profiles for the average a. 
The middle panels show the residuals of the best-Einasto fit 
for each subhalo, the thick line represent the median and 
la equivalent dispersion. The residuals are consistent with 
zero indicating that the Einasto profile fits better than NFW 
(whose residuals are shown in the bottom panel), especially 
for large a values. Interestingly, for (a) — 0.24 the NFW 
model provides a good and comparable fit to the Einasto 
model (see bottom-right panel). This is actually expected, 
since a = 0.22 represents a model that nearly follows the 
NFW profile. 

3.2 MW's dSphs constraints revisited 

|Wolf et aiT] ( |2010[ > have shown that the mass enclosed within 
the half mass radius r 1( / 2 of a dynamical system can be ro- 
bustly determined as TO1/2 = 3G~ 1 (a2)ri/2, without (pre- 
cise) knowledge of its velocity anisotropy. Here (er 2 ,) is the 
light- weighted average line-of-sight velocity dispersion of the 
system. In the case of dSph galaxies, this effectively im- 
plies a measurement of the circular velocity at ri/2, which 
therefore constrains the possible family of circular veloc- 
ity curves for a given dark matter density profile. Follow- 
ing |Boylan-Kolchin, Bullock fc Kaplinghat| ( |2011b[ ) we plot 
in Fig. [2j the 2cr constraints derived in this way for the 9 
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most luminous dSphs satellites of the MW (i.e. excluding 
the Sagittarius dwarf galaxy and the Small and Large Mag- 
ellanic Clouds). Like |Boylan-Kolchin, Bullock Sz Kapling^| 
hat (2011b I here we assume these systems are embedded in 
NFW profiles, which leads to the gray band shown in the 
figure. This band is also nearly consistent with the masses 
enclosed within 300 pc as reported by (Strigari et al.|[2008| ) 
2.5 x 10 6 < m(300 pc)/M 3.0 x 10 7 . 

The advantage of using this plot is that one can di- 
rectly compare against the results extracted from ACDM 
simulations. The filled circles in Fig. [2] show the distribu- 
tion of (r max , « max ) measured directly from our simulations 
for the satellites hosting stars. The colors indicate the pre- 
dicted luminosities and the sizes correspond to the bound 
mass fraction at present day, i.e. /bound = m(0) / m(2i n f a ii) , 
where Zi n f a n is the lowest redshift at which the progenitor 
of a given subhalo was not associated to one of the main 
Aquarius halos. Here, the values (r max , w max ) have been cor- 



rected for softening length effects following the expressions 



given by Zavala, Springel & Boylan-Kolchin ( 2010 1. As high- 



lighted in the Introduction, there are important differences 
in the location of the points from the simulations and those 
derived for the dSph satellites of the Milky Way, which may 
lead to the conclusion that there is a significant problem 
with our currently preferred cosmological model. 

However, this comparison may be need to be revisited 
since we demonstrated in the previous section that an NFW 
profile is not expected to describe well the dark matter ha- 
los of satellites in ACDM. Therefore, we have computed the 
family of (r max ,t; max ) values that are consistent with the 
measurement of v c (r 1 / 2 ) for the dSph of the Milky Way, but 
now we have considered Einasto profiles. The 2a constraints 
are shown in Fig. [2] Given the freedom we have in choosing 
the extra parameter a (this is still unconstrained by the ob- 
servations) , we have plotted 3 different bands corresponding 
to a = (0.22,0.30,0.50). The a = 0.22 (solid lines) is con- 
sistent with the NFW predictions, as expected. For larger 
values of a the constraints from observations actually over- 
lap with those found in the simulations. For a — 0.5 (dashed 
line) the match is almost perfect. 

Fig. [2] shows that there is a correlation between the 
value of a and /bound i.e. the amount of stripping a subhalo 
has experienced. Very heavily stripped objects have, on av- 
erage, higher a, and deviate the most from NFW profiles. 
The black dots in this figure correspond to the location of 
isolated dark matter halos in the same f max range as the 
satellites. This confirms that such isolated dark matter ha- 
los are well-fit by NFW profiles, and that tidal stripping is 
acting on the subhalos to change the shape of their circular 



velocity profile to the Einasto form (see also Hayashi et al 



2003). 



The Milky Way does not have many very luminous 
dwarf galaxy satellites. Brighter than Fornax (My ~ 
— 13.2), only the Sagittarius dwarf and the Magellanic 
Clouds are known, and none is included in Fig. [2] Accord- 
ing to our semi-analytic model such luminous objects would 
populate the upper-right of this plot, i.e. u max > 40 km/s 
and r max > 2.5 kpc. Therefore, although the mismatch 
between the observations (with the assumption of NFW, 
Boylan-Kolchin, Bullock & Kaplinghat |2011a l and the sim- 
ulations might be alleviated using an Einasto profile, in this 
region of the diagram, the difference may be entirely at- 
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Figure 2. Constraints for the MW's dwarf spheroidals using 
NFW (gray band) and Einasto (blue curves) profiles. Points are 
the results from the six Aquarius halos, colored according to their 
predicted luminosity and sized using the fraction of mass retained 
after infall. The black dots correspond to isolated halos in the 
simulations. 



tributed to the absence of other bright dSph in the Milky 
Way. It is the region with v max ~ 20 — 40 km/s, where we 
expect to find the subhalos hosting most of the classical 
dSph according to this plot, where the difference between 
assuming an NFW or an Einasto profile needs to be taken 
into account to bring the simulations in agreement with the 
observations. 



3.3 Effects of the host halo mass 



Springel et al. (20081 have shown that the mass function 
of dark matter halos is independent of mass, i.e. that it is 
self-similar. This implies that the number of subhalos of a 
given mass scales directly with the mass of the host (al- 
though |Gao et al.|[201 2, suggest that the slope is slightly 
larger for 10 15 Mq objects). Therefore, we can expect that, 
down to a certain scale, brighter or more massive central 
galaxies will host a larger number of satellites. This is in- 
deed shown in the top panel Fig. [3J where we have plotted 
the luminosity function of all Aquarius halos. It is clear from 
this figure that this is the case, since the three heavier of the 
Aquarius halos Aq-A-2, Aq-C-2 and Aq-D-2 have 57, 62 and 
67 satellites respectively, while the lightest, Aq-B-2 has only 
36 satellites with M v ^ —5, and hence has the shallowest 
luminosity function in the faint end. It is possible to show 
that a doubling of the mass of the host halo roughly leads 
to an increase by a factor of ~ 2 in the number of satellites 
brighter than M v — —5 (Starkcnb urg et al.|in prep I. 

We thus explore the effect of host halo mass on the 
properties of our simulated satellites, by re-scaling all ha- 
los to a common value following Hel mi, White fc Springe"!] 
( |2003[ ). Because of the scale- free nature of gravity, we may 
assume that if a halo of mass Maq is scaled to have a mass 
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Figure 3. Luminosity function for the original Aquarius sim- 
ulations (top) and once they have scaled to the mass of Aq-B-2 
(bottom). For reference we have added the observed luminosity 
function for the Milky Way as derived by |Koposov et al.| | |2008| . 
Although in the scaled version the simulations follow much more 
closely the observations, some differences remain in the number 
of satellites of a given luminosity. 
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Mmw then the subhalos's masses m should be scaled as 



~ Mmw 
m = ra—— r — = am. 
M Aq p 

Naturally the distances will become 

~ 1/3 

r — fi r, 

while the circular velocity profiles 



1/3 



(4a) 



(4b) 



(4c) 



To determine the factor [i = Mmw /Ma^, we need to spec- 
ify Mmw- As discussed in the Introduction, the value of the 
total mass of the Milky Way is quite uncertain. However, 
motivated by the remarkable match between the luminosity 
function of Aq-B-2 and that of the Milky Way ( Koposov 
et al.|2008[ ), we set M M w = M vir (Aq-B-2) = 8.2 x W^M^j 



This value is consistent with many recent studies using dif- 



ferent techniques (e.g. |Battaglia et al. 2005 2006 Smith 



|et al.||2007"l |Xue et al.||2008[ |. This implies that the value 
of n ranges from unity to 2.2 at most, which implies that 
distances and velocities in the scaled simulations will be at 
most decreased by a factor of 1.3. 

We run our semi-analytic galaxy formation model now 
for the re-scaled Aquarius simulations. The resulting lumi- 
nosity function is shown in the bottom panel of Fig.[3j where 
the new predicted magnitudes are denoted by M v . It is evi- 
dent from this figure that each halo now follows much more 



1 In this paper we denote Af v ; r = M20O1 i- e - the the mass enclosed 
in a sphere with mean density 200 times the critical value. 



Figure 4. Circular velocity profiles for scaled subhalos in three 
different luminosity bins, following the absolute magnitudes of 
the nine classical dSph of the Milky Way. The subhalos are col- 
ored according to the host halo they are associated with. This 
figure shows that the number of satellites per bin, as well as the 
velocities profiles are consistent with the measurements obtained 
for the dSph. 



closely the Milky Way's luminosity function. It is impor- 
tant to note that there is still some halo-to-halo dispersion, 
which can be attributed to the stochastic nature of the mass 
assembly of each object. That is, not all the galaxies with 
the same mass are expected to have the same number of 
satellites with the same luminosity, although some form of 
statistical equivalence should be present. 

We now study more closely the circular velocity profiles 
of the subhalos hosting satellites, since previous works have 
highlighted a discrepancy between the observations and the 



simulations ( Boylan-Kolchin, Bullock fc Kaplinghat||2011b 
|Lovell et al.||2011[ ). Fig. |4|shows the scaled circular velo city 
profiles for all the subhalos hosting satellites with luminosi- 
ties in the quoted range. As in previous sections, we have 
also included the estimates for the 9 most luminous dSphs 



of the Milky Way following Wolf et al. (20101. The first 



conclusion is that our semi-analytic model places (satellite) 
galaxies of a given luminosity in the right mass (sub)halos, 
since the amplitude of the rotation curves in all cases are 
consistent with those observed. Secondly, the number of ob- 
jects per luminosity bin is in good agreement with the num- 
ber observed, as established in Fig. [3] For example, in the 
most luminous bin (top panel) the median number of bright 
satellites per halo is 3, while for intermediate luminosities 
it is 4, and for the faintest considered here, it is 4. We em- 
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Table 1. Statistical comparison of the 9 most luminous classical 
dSphs with the simulated velocity profiles in Fig. [4] 



phasize however that the range within a given luminosity 
bin is quite broad. For example for the brightest bin, the 
scaled Aq-A-2 has just 1 satellite, while the scaled Aq-D-2 
has 7. Such large variations are not unexpected, but stresses 
that strong conclusions cannot be drawn when the number 
of objects is so small as in the case of the bright end of the 
luminosity function. 

The galaxies shown in the fainter two bins agree quite 
well with the predictions of our models. There is no sys- 
tematic mismatch, with the dSph circular velocity measure- 
ments at the half-light radius lying close to the median ve- 
locity profile of the simulated satellites. An apparent dis- 
crepancy is present in the most luminous bin —13.2 ^ M v ^ 
— 11.9 in the sense that there is a larger number of subhalos 
with circular velocities above the measured values for Fnx 
and Leo I than there is below. Nevertheless as discussed 
above, this comparison is limited by statistics and affected 
by stochastic aspects in the luminosity function. 

We quantify this by comparing the observed value of the 
circular velocity at r 1( / 2 for each dSph, with the probability 
distribution function of v c calculated at the same radius us- 
ing all the subhalos that lie in the corresponding luminosity 
bin. We compute the median v c and two percentiles of the 
distribution, namely a lower (15.9%) and an upper (84.1%) 
value, which would correspond to ±lcr in the case of a Gaus- 
sian. The results of this experiment are shown in Table [l] for 
the different satellites. Note that here we have translated the 
probability into an "equivalent" No away from the median. 
The table shows that all satellites are consistent with be- 
ing drawn from the population of subhalos hosting galaxies 
found in our simulations. 

Thus far we have focused on the nine classical dSph, and 
have excluded from the analysis the Sagittarius dwarf and 
the Magellanic Clouds. One of the questions originally posed 



by Boylan-Kolchin, Bullock & Kaplinghat (2011a I is that 



there may be a hidden population of very massive subhalos 
in the Milky Way, since the circular velocities of the classical 
dwarfs are lower than those found for the nine most massive 
(at infall) subhalos in any of the Aquarius simulations. So 
far we have shown that our model predicts the satellites 
of a given luminosity to be hosted in subhalos of the right 
mass, when comparing to the classical dSph. However we 
also need to explore what happens for systems brighter than 
Fnx, and whether we indeed expect a missing population 
from our models. As expected, the scaled Aquarius halos 
show a diverse number of systems with M v < —14, ranging 



Figure 5. Circular velocities for the subhalos present in halo 
Aq-C-2 associated to satellites with luminosities M v ^ —12. The 
symbols represent observations of the satellites of M31 in the same 
luminosity range. Open symbols represent the estimated dark 
matter contribution to v c (r) when the decomposition is available 
(see text for details). 



from 2 for Aq-C-2 to 5 for Aq-A-2. A simple comparison 
to the Milky Way satellite system would suggest that we 
cannot argue that there is a population of massive satellites 
that is missing. 

As suggested by multiple authors, the dark matter mass 
of M31 is almost as twice that of the Milky Way ( e.g. |LT| 
fc White||2008l |Kallivayalil et al.|[2009] |Guo et al.| [20101 



Watkins, Evans fc An||2010| >. This would imply that M31 
should host more, and also brighter, satellites than the Milky 
Way itself. This indeed appears to be the case, as M31 has 
9 satellites brighter than Leo I (M v — —11.9) compared to 
2 (or 5 when Sgr and the Magallenic Clouds are included) 
for the Milky Way. 

Aquarius halo Aq-C-2 has M vil = 1.77 x 10 12 M Q which 
is ~ 2.2 times larger than our candidate for the Milky Way, 
making this object a good match for M31. In Figure [5] we 
show the velocity profiles of all the nine satellites in Aq-C-2 
with M v 5j —12. We have also included measurements for 
M31's satellites with luminosities in that range. It is impor- 
tant to bear in mind that these measurements have been 
derived using a variety of methods that range from Hi rota- 
tion curve s for IC10 ( |Wilcots fc Miller|1998[ ) and M33 ( |CcF| 
belli||2003[), to 3 integral dynamical modelling for NGC147, 



NGC185 and NGC2005 DeRijckeet al. (20061. For the dSph 



Andll and AndVII the method presented by |Wolf et al.| is 
used to estimate v c (r 1 / 2 ) ( Kalirai et al.|2010 l, while for M32 
the mass is derived through Jeans modelling ( |Magorrian fc 
|Ballantyne||200l| . Many of these bright dwarf galaxies are 
not as dark matter dominated within the region populated 
by the stars as the dSph, and hence a direct comparison 
to the circular velocity of the subhalos is not quite correct. 
For example, for NGC147, NGC185 and NGC2005 the dark 
matter content is estimated to be 50%, 40% and 40%, re- 
spectively (Dc Rijck e et al. 12006] ). The open symbols in Fig. 
[5] correspond thus to the dark matter contribution to the 
circular velocity as estimated by these authors, while the 
solid points represent the total enclosed mass at the given 
radius. Clearly for M32, a very compact dwarf elliptical, the 
shown circular velocity is also an upper limit for the dark 
matter contribution. 

This comparison shows that the velocity profiles for our 
most luminous satellites in Aq-C-2 are very consistent with 



the dark halos of the MW's dSphs 7 



the observations of the dwarfs in M31 over a similar luminos- 
ity range. We therefore must conclude there is no evidence 
of a missing population of very massive of dark satellites. 



4 CONCLUSIONS 

We have used the state of the art cosmological iV-body simu- 
lations of Milky Way-like dark matter halos of the Aquarius 
project, supplemented with a semi-analytic galaxy forma- 
tion model, to study the dynamical properties of the satellite 
population in the Local Group. 

We have found that the mass profiles of the subhalos 
associated to bright satellites according to our model deviate 
from the standard NFW form (see also |Stoehr et al.||2002[ ), 
and that Einasto profiles provide much better fits. The shape 
parameter a exhibits a correlation with the amount of mass 
stripped since the time of accretion, indicating that tidal 
effects may be responsible for the changes in the profiles of 



dark matter halos once they become satellites ( Hayashi et al. 
[2003] ). 

The comparison of our models to current measurements 
of the mass enclosed within the half mass radius for the clas- 
sical dwarf spheroidals suggests that they are embedded in 
dark matter halos of fj max ~ 10 — 30 km/s with a ~ 0.2 — 0.5. 
In principle, this prediction for the values of the shape pa- 
rameter a could be tested observationally. However this re- 
quires very extensive sampling of the kinematics of stars 
near the center of the dwarf galaxies. It is also necessary 
to perform more sophisticated dynamical models, that are 
free of assumptions regarding the velocity anisotropy of the 
systems. For example Schwarzschild models of the Sculp- 
tor dSph constrain the inner logarithmic slope of the dark 
matter density profile to be d log p/d log r > —1.5 (Bred- 
dels et al.|2012 1. Better constraints could be obtained if the 
sample size were increased by a factor ~ 10. 

We have also shown that the number and internal dy- 
namics of the classical dSph in the Milky Way are consis- 
tent with the predictions of the ACDM model, if the Milky 
Way's mass is ~ 8 x 10 Mq. This value well within the 
range measured using the dynamics of stellar tracers, but 
suffers from significant uncertainties. However, it is impor- 
tant to note that this low value lowers the probability of a 
galaxy like the Milky Way to host two satellites as bright as 



the Small and Large Magellanic Clouds ( Busha et al.|2011 
Boylan-Kolchin et al.|2010 l, although such systems appear 



(20111 using 



to be rare in any case, as shown by Liu et al. 
the Sloan Digital Sky Survey. 

We have also found significant scatter in the number of 
subhalos expected to host bright satellites for the Aquarius 
halos, even when scaled to a common mass of 8 x 10 11 Mq. 
For example, the scaled Aq-A-2 has five satellites brighter 
than Fornax, while the scaled Aq-C-2 has only two (mak- 
ing it consistent with our Galaxy). Therefore, care should 
be taken to draw strong conclusions from this region of the 
luminosity function since the number of objects is small and 
heavily influenced by the host mass as well as by stochas- 
tic effects associated to particular histories. Another exam- 
ple that emphasizes this point is given by M31, which is 
nearly a factor of two more massive than the MW, and also 
has a larger number of bright satellites. Just like for the 
Milky Way, our models for the satellite population of M31 



are consistent with the observational constraints on the in- 
ternal dynamics of the brighter satellites, after taking into 
account the differences in host mass. Hence we must con- 
clude that we have found no evidence of a missing massive 
satellite problem in the Local Group. 
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